SCUPHY-TH-08001 



Extracting MSSM Masses From Heavy Higgs 
Decays to Four Leptons at the LHC 

p. Huang, N. Kerstin^, H.H. Yang 

Physics Department, Sichuan University, P.R. China 610065 

Abstract 

It is well known that finding and measuring the masses of particles in the 
Minimal Supersymmetric Standard Model (MSSM) at the Large Hadron Col- 
lider (LHC) may be possible using invariant mass distributions in exclusive 
channels containing nj jets and n; leptons. We perform this analysis for the 
{nj,ni) = (0, 4) decay of heavy Higgs bosons to neutralinos, pp — > H/A — > Xi'Xj 
= 2,3,4), which then decay to four leptons and two lightest neutralinos 
Xi via on-shell sleptons. When i = j and the sleptons are degenerate, our 
Monte Carlo study shows that the LHC will be able to measure the Higgs and 
relevant neutralino and slepton masses to roughly 30%; however, if one of these 
is already known within 5%, the other three may be found to equal or better 
accuracy. This would provide the first accurate measurement of the H/A mass 
via invariant mass distribution techniques. 
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1 Introduction 



Supersymmetry (SUSY) provides one of the most attractive candidates for physics 
beyond the Standard Model (SM) anticipated at the Large Hadron Colhder (LHC). 
The Minimal SUSY SM (MSSM), in particular, offers a rich spectrum of novel super- 
particles (sparticles) which, if relevant to LHC phenomenology, will have masses in 
the range of hundreds of GeV. Yet measuring this spectrum will be far from trivial 
since the center-of-mass (CM) energy of the LHC is available to us only by partonic 
interactions and therefore not fixed from event to event: we will not be able to tune 
the CM energy to precisely scan through mass resonances. The general consensus in 
the literature is to first conduct inclusive measurements [H El El IH IS] which require a 
suitable number of rij high energy jets plus n; isolated leptons plus missing energy (as 
carried out of the detector system by the lightest SUSY particle (LSP) in R-partity- 
conserving models), permitting a gross measurement of the sparticles' mass scale, 
and subsequently specialize to analysis of exclusive channels. Performing the latter 
in as model-independent a way as possible is a challenge in light of the tremendous 
parametric freedom in the MSSM, where sparticle decay topologies vary widely across 
this parameter space. 

Once a specific decay chain can be identified we can construct relativistically- 
invariant combinations of rij and rii momenta and analyze their distributions over 
a large number of events. Since endpoints of these distributions are typically well- 
defined analytic functions of sparticle masses for this given decay topology, a suf- 
ficiently large and pure sample of sparticle events can provide a clean endpoint 
which constrains some set of sparticle masses. Examples of situations where this 
technique yields promising results include where we can identify X2 ^^^^Xi {one 

endpoint for two unknown sparticle masses), II — ^ X2^^ ~^ ^'r ^'^^^ ~^ I'^l'^l'^Xi 
(one endpoint for three masses), pp — > Ql^l with qi qX2 ~^ l^l^q — > l^l^qxi 
or qi — » qx2 hqxi (six endpoints for four masses ), and g qq ^ X2II ~^ 
l^r^qq — i> l^r^qqxi (seven endpoints for five masses), all of which are discussed in 
more detail in Refs. [H], [Z|, and [S]. However, in these examples either the number 
of endpoints is less than the number of unknown masses, so we can only constrain 
the MSSM to some surface in mass parameter space, and/or we have to make extra 
model assumptions which ensure the assumed decay topology occurs with a sufficient 
rate. 

In this paper we investigate whether this technique performs better with a larger 
number rii of final state leptons. As signals consisting of isolated leptons provide 
cleaner signals than those from jets at a hadron collider, invariant masses constructed 
from the former may be subject to less error. While up to now researchers have con- 
sidered rz; = 3 which only gives three independent invariant combinations, higher 
values of ni rapidly give more: n; = 4 already gives seven combination^. For virtu- 
ally any conceivable MSSM decay chain this would already suffice to overconstrain 

^For an endstate with rij jets and ni leptons {N = nj + n;), the number of independent invariant 
combinations is equal to the number of pairwise contractions, 'N choose 2', plus the number of 
contractions with e^j^po-, 'N choose 4'. iV = 4, 5, and 6 gives 7, 15, and 30 combinations, respectively. 



the unknown masses if the endpoints of these distributions were precisely known. 
However, with this rise of the number of constraints comes the lower statistics on 
each constraint since cross- sect ions naturally fall with rii. But backgrounds (both 
SUSY and SM) will correspondingly become smaller, hence confidence in the as- 
sumed decay channel will increase. Moreover, there is information in the peaks of 
the distributions as well - by definition these have higher statistics and therefore 
lower error, though we would have to investigate the sensitivity of these (or any 
other 'shape' variable) to cuts. A complex interplay of all these factors therefore 
determines the success of this programme, which ultimately depends on the specific 
endstate considered. 

In this work we concentrate on the {nj = 0, ni = 4) endstate which may result 
from the decays of heavy Higgs bosons and (hereafter collectively referred to 
as 'Higgs') to neutralinos, 

PP - H/A ^ - k^'h^k^k^ ^ /i^/'i'=/2^/'2^X?X? (1) 

proceeding via on-shell sleptons of electron or muon (/i^2 C {e, fi}) flavor. As shown in 
Ref. [H] this signal has favorable rates (> 100 events for 100 fb~^ integrated luminosity 
at the LHC) over much of the {fi, M2)-plan^ when the Higgs mass is in a favorable 
range (350 GeV < Ma < 700 GeV) and flrst and second generation slepton masses are 
sufficiently light (< 200 GeV). Here the correlation between lepton invariant mass 
pairs, i.e. a wedgebox plot [IHl [H], tells us something about the mass differences 
— and m^o — m^o: they are zero for a symmetric boxlike wedgebox plot. In 
this {i = j) case the only significant background after a suitable jet cut consists of 



charginoci, i.e. pp xfxf j which may subsequently decay as 
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Among such four-lepton endstates, the ratio of flavor-balanced (e+e~e"'"e~-|-/i"'"/i~/i"'"/i~ + 
e+e~/i+/i~ ) to flavor-unbalanced (e"'"e~e^/i'^ + ix'^ ix~ im^ e"^) events, which we deflne 
as 7^-1- , is close to unity since there is no correlation between the flavors of I and 
I" . On the other hand, a pure Higgs signal would have no flavor-unbalanced events 
and hence TZ± oo. A boxlike wedgebox plot, therefore, not only contains much 
information in leptonic invariant mass distributions (potentially seven upper and 
seven lower endpoints, if these are indeed independent, plus constraints from seven 
peaks) which can be used to flnd the six unknown masses with high precision, but 
also provides a model-independent confldence measure {JZ±). 

The rest of this paper is organized as follows: in Section [2] we present our deriva- 
tion of analytical expressions for the endpoints of the seven invariant combinations 



■'Gauge unification here fixes Mi = 5/3 tan^ OwM2. 

*One might worry about neutralino channels pp ^ X^X% however, as shown in Ref. |10j . these 
require a couphng to the Z-boson which, in the notation of [T5], is proportional to {Ni^N*^ — NiiN*^); 
here the crucial minus sign arises from the different hypercharges of the two MSSM Higgs doublets. 
For i = j this coupling is highly suppressed; i ^ j processes can be significant, but these do not give 
a boxlike wedgebox plot. 



of lepton momenta (since exact formulae are rather lengthy they are collected in the 
Appendix); in Section [3] we test this method with Monte Carlo (MC) simulation of 
LHC data generated at three different MSSM parameter points with backgrounds 
and detector effects included. Section H] summarizes and discusses these results. 



2 Endpoint Theory 

From the four-lepton endstate of ([1]) depicted in Fig. [H we can form six independent 
relativistically invariant bi-contractions (we take lepton masses to be zero in the 
following) , 

{piripv),, {PinP2),, {PiY{P2'),. {pvY{P2),. {Pi'nP2'),, iP2nP2'),i2) 

in addition to the totally antisymmetric invariant 

^4 = PlPl'P2P2'(^liupa (3) 

However, when we construct a distribution of a function of these invariants we must 
take care that this function is totally symmetric under interchanges of labels 1^1' 
and 2^2' because of ambiguity in lepton identification. The usual dilepton invariant 
masses 

M21 = 2{pinpy)^ , M2, = 2(p2)^(P2')m (4) 

obey this rule. 

But because we want to include same-flavor endstates {e'^e~e~^e~ and n'^ ii~ [i^ 
in addition to those of opposite flavor (e"'"e~/i"'"/i~), we choose a list of positive def- 
inite functions of the invariants in ([2]) which are symmetric under all label inter- 
changes: 

Ml = {pi+Pv +P2+P2'f (5) 
= {{Pl+PV - P2- P2'Y + {Pl+ P2' - P2- Pl'Y + {P1+P2- PV - P2'Y}/^ 



= {{Pl+PV +P2-P2'Y + (Pl+PV +P2' -P2, 

+ (Pl +P2+ P2' - Pl'Y + {P2 + P2' + PV - Pl)^}/4 

= {{Pl + PV ~ P2Y + {Pl + PV - P2'Y + ij>l + P2 - P2'Y + {P2 + P2' - Pl'Y 

+ {Pl - PV + P2Y + {Pl - Pi' + P2'Y + (Pl -P2+ P2'Y + {P2 - P2' + PvY 

+ {PV +P2- PlY + {PV + P2' - PlY + {P2 + P2' - PlY + iP2' + Pl' - P2)^}/12 

= {(p, +p2)^ + (pi+Pl'+P2')^+(Pl+P2+P2')' + (P2+P2'+Pl')'}/4 

M^l = {(pi + prY + (Pl + P2'Y + (Pl + P2Y + iP2 + P2'Y + iP2 + Pl'Y + iPl' + P2'Y}/Q 



^Symmetric linear combinations of ^ can only be proportional to the four-lepton invariant mass 
M41. To get other independent combinations we must, in the simplest scheme, choose averages of 
squares of dl]). In doing so all lower endpoints become zero. 



Figure 1: Kinematic picture of the Higgs decay showing angular variables defined in the 
rest frames of the decaying particle at the corresponding vertex. 

With this set of invariants there will be no combinatoric background from lepton 
misidentification; on the other hand, since the minimum values of ([5]) are in fact 
zero, we only have upper endpoints (hereafter simply 'endpoints') to measure; the 
former effect should be more of an advantage than the latter a disadvantage. Since 
the seven independent endpoints of and are well-defined functions of the six 
masses rriA, frif^, m^-, ""^jjo, ttl^o, and fn^, endpoint measurements would plausibly 
lead to a system of equations which overconstrains these masses. 

To derive analytical formulae for the endpoints, we start by defining the kinematic 
degrees of freedom of the decay chain ([1]) in Fig. [TJ Since all decays are two-body, 
daughter particles are produced back-to-back in the rest frame of the decaying par- 
ticle with spherical angles {6, 0), distributed evenly in the ranges < (p < 2n and 
—1 < cos 6' < 1 (we ignore spin effects; this will be justified by comparison with 
Monte Carlo results below). The momenta and directions of all four leptons and the 
LSPs are then completely determined by the set of 8 angles {9i, {9j, (pj), {6i, 0i), 
and (6*2, 02) which are defined in the rest frames of the /i, and I2, respectively 

(for exact procedure, please see Appendix). 

Since (E]) define a set of invariant masses, they will in general be maximized 
when some particular sum of lepton energies is maximized while a corresponding 
sum of 3-momentum is minimized; this type of situation is realized when each 
polar angle 6i, 6j, 6'i,and 62 equals or vr, values of the various angles then 
being irrelevant. Each expression in therefore has potentially 2^ = 16 possi- 
ble endpoints (though some of these are the same by symmetry), the maximum 
(true endpoint) of which depends on the precise values of ratios among the masses 
niA, rrif^, rrif^, m^o, m jjro , and m^o . 

To illustrate this phenomenon, consider M41 with i = j = 2 and I2 = h = I for 
simplicity. Depending on the hierarchy of the ratios 



r2A = 




Table 1: Correct endpoint expressions for the four-lepton invariant mass M4; for all possible 
orderings of the ratios defined in 
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we find that one of the following three expressions is the endpoint: 

M,i[--++] ^ , (7) 



M.,\+ + 1 = 
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where the notation suggests the values of cos^j, cos 6j, cos^i, and cos02 which 
equal 1 ('+') or -1 ('-'). Each of the six possible orderings of the ratios in ([6]) 
single out the largest of these endpoints, as summarized in Table [H Thus, kinematic 
configurations that maximize M4/ change as some masses become larger relative to 
others. In the particular limit where niA is much larger than the other masses, the 
hierarchy listed in the first row of Table [1] applies with the angular configuration of 
Fig. [2^; this makes physical sense since the leptons should like to maximize their 
energies by being emitted parallel to their highly energetic mother neutralinos. As 
the Higgs mass is lowered, however, it may maximize total leptonic energy to emit 
one or both leptons antiparallel (cf. Fig. [2)d,c) depending on the precise ratios (jS]). 

A similar situation applies to the other invariants in ([5]): there exist three different 
endpoints depending on mass ratios, though the physical interpretation is not as 
clear as for Mai. We have collected analytical expressions for all such endpoints in 
the AppendixO 

The distributions of ([3]) and ([5]) also contain information in their shapes and 
specifically their peaks. We have not derived analytical formulae for these (though 
in principle this is possible via the method of Ref. [13]). Peak values of distributions 
can always be numerically computed and we will find these useful as additional 



^We have not derived an endpoint for 04 because the tail of this distribution is extremely shallow 
and therefore unlikely to be useful in the MC analysis. 
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Figure 2: Angular configuration for maximizing the four-lepton invariant mass in the 
limit where the Higgs mass is very large (a), or for lower Higgs ' masses with < r^s 
(b) or ru > (c). 

constraints (this is in fact the only use we have found for a^). A more detailed 
discussion of this strategy appears in the following sections. 



3 Monte Carlo Results 

In this section we would like to investigate how well the programme sketched above 
works for Monte Carlo (MC) LHC events simulated at actual points in MSSM pa- 
rameter space. We employ the HERWIG 6.5 [14] MC package (whose MSSM input 
information comes from ISASUSY [I5] through ISAWIG [I6] and HDECAY [n]fl 
to generate LHC events for an integrated luminosity of 300 fb~^, roughly equivalent 
to several years' high luminosity data, which we then run through private programs 
simulating a typical LHC detector environment. 

We pass pp H/A events with four hard and isolateci§ leptons with flavor 
structures e"*"e~/i+yU~, e'^e~e'^e~, or fi'^ fi~ fi'^ fi~ (hereafter designated the 'isolated- 
41 cut') at the following three parameter points: 

^ The CTEQ 6M [18] set of parton distribution functions is used with top and bottom quark 
masses set to nit = 175 GeV and mf, — 4.25 GeV, respectively. 

®SpecificaUy pfj, > 10, 8 GeV for e"*",^^, respectively; \ri^\ < 2.4; isolation requires no tracks of 
other charged particles in a r = 0.3 rad cone around the lepton, and less than 3 GeV of energy 
deposited into the electromagnetic calorimeter for 0.05 rad < r < 0.3 rad around the lepton. 



Table 2: Relevant physical masses at the parameter points (all masses in GeV). 









xll 


X^ 


ei, /ii 


AO 


i/o 


MSSM 1 


127.7 


243.9 


415.3 


436.1 


156.2 


600.0 


600.4 


MSSM 2 


120.9 


158.1 


205.4 


399.9 


156.2 


500.0 


500.4 


SPSla 


96.1 


176.8 


358.8 


377.8 


143.0 


393.6 


394.0 



MSSM Sample Point 1 

/i = 410 GeV Ms = 260 GeV tan /3 = 20 
Ms4iL,R = 150 GeV M^^^ = 250 GeV 
Ma = 600 GeV Mg^g = 900 GeV 

MSSM Sample Point 2 

H = 150 GeV M2 = 380 GeV tan p = 10 
Ms;-^^ ^ = 150 GeV M^^ ^ = 250 GeV 
Ma = 500 GeV Mg;^ = 900 GeV 

Snowmass Benchmark SPSla |19| 

/i = 352 GeV M2 = 193 GeV tan (3 = 10 
Mg;-^^ « = 136 GeV M^^ ^ = 135 GeV 
Ma = 394 GeV ^ 500 GeV Mg ^ 600 GeV 

with other relevant sparticle masses at these points appearing in Table [2l 

These points are chosen because the dominant source of signal events is H/A — > 
3^o^o _ 2 for MSSMl and SPSla, i = 3 for MSSM2). This can in fact be deduced 
from kinematics alone as follows: in Ref. [10] it was demonstrated that all box- 
like wedgebox plots with low jet activity follow predominantly from either Higgs or 
chargino pair decays; the dominating presence of the former over the latter can be es- 
timated by the ratio of flavor-balanced to flavor- unbalanced events, TZ±, as explained 
in the Introduction. Since all three of these points have a 'simple box' wedgebox plot 
with TZ± » 1, we could know prior to knowing any MSSM parameters that the 
observed four-lepton events were mostly from the Higgs decays ([T]) with i = j and 
degenerate slepton masse^. 

Before conducting a full investigation with backgrounds and further cuts, we first 
investigate the signal pp —>■ H/A channel only at Sample Point 1 so as to understand 
how well endpoints and peaks agree with their theoretical values. 

^We could also deduce that the sleptons are on-shell since off-shell decays have very different 
invariant mass distributions. See Appendix for discussion. 




Figure 3: Monte Carlo results for H/A signal only at MSSM Sample Point 1 for 300/6^ , 
with spline fits to theoretical distributions (Id^ events from kinematics only) overlaid and 
normalized (by signal rate/10^). 

3.1 H/A Signal Only 

The dominant source of lepton pairs at Sample Point 1 is 

pp ^ H/A ^ ^,xl ^ h^h^h^k'^ - /i^/;^/2^r2^X?X? (8) 

For 300 fb~^ luminosity we obtain nearly 1000 H/A events after the isolated-41 cut; 
this rate is quite good but by no means optimized. With a lighter Higgs or higher 
tan/3 (increasing the Higgs production cross section) or heavier L-handed sleptons 
(reducing spoiler modes such as X2 ~^ which compete with decays to e,jl), for 
example, the rate could increase by a factor of several [TO]. Distributions of all 
kinematic invariants are shown in Fig. [31 We first notice that the wedgebox plot 
is a very crisp symmetrical box bounded by the M21 edge ~ 108 GeV as expected. 
Other plots have overlying normalized fits to purely theoretical distributions which 
we have computed from relativistic kinematics only. We see a reasonable agreement 
between the two which, of course, is not expected to be perfect since the theoretical 
distributions do not include the isolated-41 cut, detector effects, particle widths, and 



Table 3: Comparison of Endpoints (EP) and Peaks (P) of Monte Carlo (MC) and theo- 
retical values (boldface, based on relativistic kinematics only), for signal only at MSSM 
Sample Point 1 (all numbers to GeV precision). Error is quadrature sum of bin size, fit 
interval, and fit parameter error. Error on theory peaks arises from using a Gaussian to fit 
a slightly non-Gaussian shape. (*)The endpoint for \ai\^/'^ is very shallow and hard to fit, 
so the corresponding entry is blank. (**) We do not fit the peak for M21 which theoretically 
coincides with its endpoint. 





M21 


Mil 


M2121 




M121 




Mu 




EP 


105 ±5 


340 ± 10 


275 ± 10 


285 ± 10 


230 ± 10 


262 ± 10 


130 ± 10 


* 




108 


345 


282 


288 


235 


266 


131 


* 


P 


** 


228 ± 17 


142 ± 10 


131 ± 15 


116 ± 10 


167 ± 11 


72 ±5 


63 ±10 




** 


230 ±2 


142 ±2 


132 ± 5 


108 ±3 


168 ±2 


72 ± 1 


60 ±1 



other features of the MC. Presumably such account for the suppressed peak of M21 
and elongated tail of |a4|^/^, neither of which we will employ in any fits for this work. 

Peaks and tails of MC distributions are generally well-fit with Gaussian func- 
tion£^. However, theoretical peaks may occasionally (as in the case of M^a; in Fig. [3]) 
be somewhat asymmetric which gives rise to a slight (2 — 5 GeV) error in compari- 
son to the MC fits, deemed tolerable at this level of analysis. We fit all histograms 
interactively in PAW, identifying three sources of error: binning, fit interval, and 
statistical error on the parameters of the fit; these are added in quadrature. The 
small error in fitting theory peaks then arises almost totally from fit interval error. 
For MC data, endpoints are dominated by binning and fit interval, and peaks by fit 
parameter error. The agreement between fits for the signal and theory is very good 
(cf. Table [3]) and stable when subject to a liberal range of missing energy and jet cuts 
(Table Hj), the use of which we will find necessary when backgrounds are included. 

3.2 With Backgrounds and Cuts 

Confident that a pure signal can be fit closely to the theory with a simple fitting 
procedure, we now add all SM and MSSM backgrounds into our analysis. SM back- 
grounds (tt, th, tiZ^, W~^W~, Z^Z^) are eliminate(£ll by the isolated-41 cut in addition 

Specifically, we fit peaks to a 3-parameter function a + 6 exp[—{x — c)^] (a parabola would work 
just as well). Although tails happen to look quite linear at Sample Point 1, we have found that this 
is not generally true. There is always a slight amount of nonlinearity which, for economy's sake, we 
fit to the same functional form as used for the peaks. 

Residual pp z(*)z(*) (i.e. either Z-boson may be on- or off-shell) backgrounds can be large 
(up to 1000 events after the isolated-41 cut) but are constant in MSSM space. In our analysis 
we generated six backgrounds, using one subtracted from the average of the other five; 

cancellation was nearly total. Where this background, more-or-less forming two strips of events 
within T z of mz on the wedgebox plot, is external to the box structure, it can also be subtracted 
by inspection. 



Table 4: Effects of missing energy (I^) and jet (Ej) cuts on fitted values of peaks and 
endpoints (P\EP) for the signal only at MSSM Point 1 (all units are GeV). Variation of 
peak and endpoint values are seen to be about ±5 GeV and ±3 GeV, respectively. This is 
smaller than the fitting error on these. 





Mil 


M2121 
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Mu 


04 


1/4 


(0,0) 


241 


345 


150 


272 


124 


275 


108 


227 


170 


264 


72 


128 


66 




(0,30) 


244 


347 


152 


272 


125 


274 


114 


227 


173 


264 


74 


128 


64 




(0,50) 


240 


344 


151 


273 


134 


277 


117 


228 


169 


262 


76 


129 


68 




(0,100) 


241 


342 


152 


278 


137 


276 


115 


232 


173 


265 


75 


128 


69 




(10,0) 


238 


345 


154 


271 


132 


275 


119 


228 


175 


262 


74 


131 


68 




(10,50) 


242 


343 


148 


278 


132 


275 


117 


227 


174 


264 


75 


128 


66 




(10,100) 


239 


345 


152 


272 


138 


280 


118 


236 


177 


266 


74 


127 


65 




(30,0) 


231 


344 


146 


273 


136 


281 


115 


229 


171 


264 


74 


128 


68 




(30,50) 


232 


346 


149 


270 


136 


278 


115 


231 


166 


260 


75 


129 


61 




(30, 100) 


233 


346 


150 


271 


129 


282 


116 


230 


175 


264 


74 


130 


66 





to demanding at least > 20 GeV missing energy in each event. 

The main SUSY backgrounds at Sample Points 1 and 2 are gaugino processes 
{(lj9)Xi'^: x'iXp and xfxj: plus a handful of events from tih, and tH~ (the 
reader may refer to our earlier works [lOl [H] for details on all backgrounds and 
cut efficiencies for this four-lepton signal). Backgrounds with colored sparticles can 
be suppressed via a jet cut: only events with jet^ having reconstructed energy 
Ej < 50 GeV are allowed to pass. Chargino and slepton backgrounds are greatly 
reduced by subtracting distributions of flavor-unbalanced events. 

In Fig. IHwe see the resulting distributions at Sample Point 1, which after cuts 
has a large value of 7l± = 786/48 ~ 16.3. The wedgebox plot (Fig. [5]) retains a very 
dense box shape, though now with a few dozen background events scattered outside 
the perimeter. Peaks and endpoints have drifted but they are still in good agreement 
with theoretical values (see TableOfor results at this and the other parameter points). 
The drift is mostly due to backgrounds, since the missing energy and jet cuts only 
have a mild (±5 GeV or less) effect when these are varied across a liberal range (cf. 
Table m). 

At Sample Point 2 the same backgrounds are more formidable, though here TZ± = 
197/26 ~ 7.6 is still large and we have a clear boxlike wedgebox topology in Fig. O 
Therefore, by our hypothesis fitted values should be in reasonable agreement with 
those for a Higgs decay H/A y°Y° (note, however, that this technique does not 
tell us which neutralino is involveq^j). As the shapes of the various distributions look 

"'^^Jets are defined by a cone algorithm with r — 0.4 and must have \ri^\ < 2.4 . 

^■^Here, in fact, i = 3 since the smallness of the HX2X2 coupUng is more severe than phase-space 
suppression to the heavier Xg-pair. Illustrating this interplay of parameters forms another reason 
for our choice of this MSSM parameter point. 



Table 5: Comparison of MC Endpoints (EP) and Peaks (P) (all numbers to GeV preci- 
sion)for MSSM Sample Points 1,2 and SPSla to theoretical values in boldface. Theoretical 
endpoints are calculated from formulae in the Appendix, whereas all other values are fits as 
described in the text. 





MSSM1 -FP 


MSSM1 -P 




MSSM9-P 
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107 ± 5 


** 


83 ± 5 


** 


75 ± 5 


** 




1 08 


** 


85 


** 


77 


** 


Ma, 


345 ± 10 


235 ± 23 


262 ± 10 


182 ± 16 


212 ± 20 


143 ± 17 




345 


230 ±2 


257 


183 ±2 


200 


148 ±2 


M2121 


275 ± 10 


141 ± 10 


185 ±15 


109 ± 10 


140 ± 20 


86± 11 




282 


142 ± 2 


196 


112 ± 2 


156 


91 ± 2 


Mm 


290 ± 10 


131 ± 16 


163 ±15 


89 ± 11 


160 ± 20 


69 ±10 




288 


132 ± 5 


180 


89 ±5 


151 


72 ± 2 


M121 


240 ± 10 


111 ± 10 


139 ±15 


86 ±10 


135 ±20 


66 ± 10 




235 


108 ±3 


147 


89 ±2 


129 


71 ± 2 


M^i 


271 ± 10 


164 ± 11 


190 ± 10 


132 ± 11 


150 ± 20 


103 ± 10 




266 


168 ±2 


186 


134 ± 2 


148 


107 ± 2 


Ma 


135 ±6 


72 ±5 


86 ±8 


57 ±5 


90 ±20 


44± 10 




131 


72 ± 1 


86 


58 ±1 


72 


46 ± 1 




** 


60 ±8 


** 


50 ± 12 


** 


41 ± 10 




** 


60 ± 1 


** 


42 ± 1 


** 


35 ±1 



rather similar to those at Sample Point 1, we do not display them. Lower statistics 
lead to larger fitting errors, though agreement with theory is nevertheless intact. 

Snowmass Benchmark Point SPSla would seem to present a more challenging 
case. With a lighter spectrum of colored sparticles, there arise significant back- 
grounds due to squarks and gluinos which force us to tighten the jet cut to Ej < 30 GeV; 
this however diminishes the signal {H/A X2X2) close to the level of the irreducible 
SUSY backgrounds, giving a lower value of TZ± = 138/59 ~ 2.3. Statistics now barely 
allow for the identification of a box topology in the wedgebox plot (Fig. E]) and we 
would seem to be testing the limits of our technique in fitting the distributions in 
Fig. [6l Yet fitted central values are reasonably close to theoretical expectations in 
Table [5] and still agree within error-bars. 

3.3 Mass Extraction 

We now wish to extract the masses rriA, ttlj, m^o, and m^o and their associated 
uncertainties from the MC data in Table [51 First we search in mass space for solutions 
(within errors) to the seven endpoint constraints using the formulae in the Appendix; 
in so doing we assume that the M21 edge has negligible error since if there is any 
modest rate for X2 ~^ Xi^^^"^ via Higgs channels then surely there is a much greater 





Figure 4: Monte Carlo results (for 300 fb~ ) at MSSM Sample Point 1 for signal and all 
backgrounds with missing energy and jet energy cuts described in the text. Flavor unbalanced 
distributions have been subtracted in this and the following plots. 

rate via squark^, e.g. qi l^V^q l^l^q^-, so this edge will have been 

already measured very precisely (perhaps at the 0.1 GeV level [SO])- We thus use the 
M21 formula to solve for m^o while scanning over m^, mj, and m^o for combinations 
that satisfy the other six endpoints within their assigned errors from Table [51 As 
the reader can verify from the Appendix, the endpoint functions compose a highly 
nonlinear system of equations which, when each endpoint value has a small (percent 
level) uncertainty, generally gives rise to a discrete set of solutions (perhaps eight). In 



We may infer this, for example, from the presence of very jetty backgrounds. 
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Figure 5: Wedgebox Plots (for 300 at MSSM Sample Point 1 (left), Sample Point 2 
(center), and SPSla (right) for signal and all backgrounds after cuts. 

our case the uncertainty on each endpoint is somewhat larger, causing these discrete 
solutions to merge into a continuous range as shown in the first row of Table El 
However, many solutions with the same endpoints have very different (by say, tens 
of GeV) peakO — though our peak resolutions are not always able to resolve these 
at SPSla and Sample Point 2, we may narrow down the range of solutions somewhat 
at Sample Point 1 as seen in the second row of Table [61 If, in addition, the LSP mass 
is already known to fair precision {e.g. ±5 GeV) from inclusive measurements [21] 
or other methods then all masses at Sample Points 1 and 2 can be constrained well 
within 5% of their nominal values. SPSla does not fare as well (except for the Higgs 
mass, which is also within 5%) but there are other techniques to constrain masses 
at this well-studied MSSM point which, when combined with ours, would certainly 
yield higher precision. 

4 Summary and Discussion 

In the foregoing we have investigated seven invariant mass distributions {M41, M21211 
Mi3i, M121, M31, Mil, and 04) of the four-lepton endstate of the Higgs' decay ([I]), 
finding analytical expressions for endpoints of six of these {M41, M2121, Misi, M121, 
M^i, and Mu) in the case of degenerate neutralinos {i = j) and sleptons. Each 
invariant has three possible endpoint expressions depending on the precise ratios of 
sparticle masses; though we explicitly found these ratios for M41, in practice one 
simply uses the maximum of the three endpoints as the correct one. Adding the 
well known dilepton invariant mass edge M21 to this list then gives seven invariant 
functions of rriA, fnj, m^o, and fn^o. We note this list of invariants is not unique: 
one could in principle select a totally different set (though still functions of ([2])) 
which might give sharper endpoints or peaks and improve upon the results we have 

^^We know that theoretical peaks are always close to MC peaks (compared to error-bars) so this 
method will suffice; in a comparison between actual LHC data and theory we would use a more 
sophisticated estimate of peaks (at the level of the best available MC), or fit to whole distribution 
shapes. 





Figure 6: Monte Carlo results (for 300/6 ) at SPSla for signal and all backgrounds with 
cuts described in the text. 



presented here. 

One advantage of concentrating on the Higgs decay ([T]) is that one need not 
assume any specific choice of MSSM parameters: the signal alone, i.e. a boxlike 
wedgebox plot with a high ratio Tl± of fiavor balanced/unbalanced events identifies 
the responsible decay^. One might then expect that any data sample of Higgs de- 
cays ([1]) with degenerate neutralinos and sleptons giving even modest determination 
of the seven endpoints would allow one to overconstrain the four sparticle masses 



^^Conversely, if the signal is absent then we can rule out a significant region of MSSM parameter 
space. 



Table 6: Extracted values of masses at MSSM Points 1,2 and SPSla from matching theory 
to endpoints only ("EP"), endpoints and peaks ("EP+P"), or these in addition to constrain- 
ing the LSP mass to within 5 GeV of its true value ("C.LSP"). 
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EP 


mj 
rriA 


136 ± 54 
171 ±80 
284 ± 88 
627 ± 123 


144 ± 46 
185 ±51 
230 ± 47 
551 ± 102 


135 ± 55 
177 ± 88 
227 ± 70 
473 ± 121 


EP + P 


Al 

mj 

Ai 

niA 


137 ±53 
164 ± 56 
255 ± 56 
623 ± 119 


144 ± 46 
185 ±51 
230 ± 47 
551 ± 102 


135 ± 55 
177 ± 88 
222 ± 65 
473 ± 121 


C.LSP 


Al 

rriA 


128 ±5 

155 ±7 
242 ±7 
603 ± 16 


121 ± 5 

156 ±5 
208 ±3 
497 ± 11 


96 ±5 

146 ± 37 
188 ± 19 
393 ± 17 



rriA, TTLi, m^o, and m^o. Were this determination more precise (at the GeV level) 
this would be so, at least up to a n-fold discrete ambiguity, but in realistic MC sim- 
ulation we find these endpoints typically have a precision on the order of 10 GeV or 
more — in this case masses can be constrained only within a large (50 — 100 GeV) 
continuous range. Results improve with inclusion of peak constraints, but accuracies 
still linger at the 30% level on average. If one of the masses {e.g. ttl^o) is already 
known with some precision then some or all of the other masses can be found with 
equal or better precision, depending on statistics and purity of sample (7?.±). In this 
way. Sample Points 1 and 2 lie in the center of disjoint regions in the (/i, M2)-plane 
where percent-level determination of all the masses is possible. More sophisticated 
analysis fitting to whole distribution shapes could in principle be employed for supe- 
rior results. We could also couple this method with other techniques in the literature 
which do not employ invariant distributions [221 [23]; such a 'hybrid' [21] analysis, 
particularly suitable to SPSla for example, would undoubtedly give much better 
mass determination. 

Our analysis could also be extended to colored sparticle decays where rates and 
therefore statistics are much higher. Here we have more avenues to explore: the 
gluino and squarks typically decay through multiple competing channels via more 
intermediate states, demanding a more intricate analysis; but we have jets in addition 
to leptons, so it may be possible to consider higher values of = rij ± ni, giving 
many more invariants. The {nj,ni) = (2,2) work of Ref. [BJ is a promising step in 
this direction. In this connection, the concept of the wedgebox plot, i.e. a plot of di- 
electron versus di-muon invariant masses, could be extended to correlations between 
any two or three invariants which might prove as useful in isolating specific colored 
sparticle decays as we have found it to be invaluable in identifying the Higgs decays 



under study in the present work. 

Though MSSM parameters will ultimately be determined via a global fit to all 
available data, it is nevertheless important to have some model-independent idea 
of the rough values of these parameters as a starting point. We envision a general 
strategy of considering all possible N = rij + rii final states' kinematic invariant 
distributions — endpoints, shapes, and correlations among these — to not only 
single out specific decay chains but also pinpoint the masses in these. While up to 
now the community of SUSY phenomenologists had fruitfully executed this strategy 
up iV = 4 (n; = 2), in this work we have pushed the frontier to = 4 (n; = 4) with 
encouraging results. 
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Appendix 



We outline here the derivation of endpoint expressions (with the exception of 04); 
to do this, we need to compute the four-momentum of each of the four leptons as 
functions of the eight spherical angles (defined in Fig. [T]) in the frame of the decaying 
Higgs, and then apply ([5]). 

Let us start with the selectron, denoted h in Fig. [H In its rest frame, the LSP 
and electron are produced back-to-back: 



Pi 



( E, \ 
Pi sin 9i cos 01 
Pi sin 61 sin (pi 
Pi cos 6*1 



p: 



( E^- \ 

—pi sin 6*1 cos 01 
—pi sin 6*1 sin 0i 
—pi cos 61 



where we have imposed momentum conservation 
El = i , =m,-Ei 



2m„ 



~Pi + ~P 

\Pi\ = 



0) and we also have 



2 2 
— ml 

2mo 



in an abbreviated notation: m^ = m^o , mi = wi^^o, and = m^^. Now we back- 
wards Lorentz boost the electron along the z-direction by the selectron momentum 
(defined in the rest frame of the decaying x?); the parameters of this boost are 
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and in this frame the positron has momentum 
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We now have the positron and electron momenta in the frame of the parent neu- 
tralino. To get these in the frame of the Higgs, we first rotate the coordinates so that 
the neutralino's momentum is along a new z-axis; this requires a rotation matrix of 
the form 

cos 9i cos (pi — sin (pi sin 9i cos ( 
cos 6i sin (pi cos (pi sin 6i sin < 
— sin 9i cos 9i 

And now we may backwards Lorentz boost by the neutralino momentum along z 
using 

mj — m^ + m\ 

li = 7, ) o 

zmimA ZmimA 

The muon and antimuon momenta can be found analogously (being attentive to sign 
changes, of course) with the obvious substitutions of i — >• j, 1 — *■ 2, and m^ — >• mg'. 



^{mj - m^ + m\Y 



4mfm\ 



Endpoint expressions (jl]) and are maximal when polar angles cos^j, cos6j, 
cos^i, and cos ^2 are equal to ±1. For the well-known dilepton mass edge, we of 
course have 

Other endpoints have three maxima (labelled by polar angles 1 ('+') or -1 ('-') in 
brackets [±±±±] as described in the text) which, for the case we consider where i = j 
and sleptons are degenerate, are functions of niA, rrii, rus, and mi. The following 
definitions assist in writing each expression more compactly: 



A = \/m\ — 4:m\mi 



J = \lmf — ml 



= - i rnA ± \/m\ — Amf 



Hu = a ml -b ml 



M4i[+ ] = i^~r^^ {mlmjiml + A) - A + m\{-2mjml + mfj) ; 

M4i[ ] = 2mAU^^ [mlmj{-m\ + A) - A + m\{2m1ml - )) ; 

M2i2i[- + +-] = ^^^^^{3mf + m^^Hl^ + 16m]miHf + 32m8 + 4mfHl + 6mf{m]Hl + 
8mf) + 4mj{mfHl^ + mfif24)|i/4. 

M2i2i[+ H ] = J/4 i 2 {3mfm^F+-4mfmfm^(3F+-J^[A-2mf+m^])-4m^m|mg(3m1.m|+ 

32mf + 24m^mf - m^A - 8mf [A + 6ml] + [3mf A - 2mf] + m^[mf - 16mf - Smfmj + 
3mf (A + 2ml)] ) + Gmfmfmj (3F+ + [lOmf - Gm^ - 2 A - 2m\] ) + {3m\ \m\ - A] + 
3m\mf[m\ + A] - m^[32mf + 52mfmp'^ + 12mf m^] + 2mf [32m,^° - 64mfml + 3mf A - 
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32mf + 24mfml + A + 8mf [A - Gm^] - m2[3mf A + 2mf] + m^[m^ - 16mf - 3mf - 
3mf (A - 2ml)]) - 6mf mfmf ( J^flOmf - Qml -2A- 2m\] - 3F-) + m^,{3m\J^[m\ + A] + 
3m\mf[m\ - A] - m\[32m^J'^ + 52mf + 12mfm^] + 2mf [32mf - 64mfml - 3mf A + 
16mfml{A + 3ml) - 8mf (3mf A + 2mf ) + mf{3ml + Mm^A)])}^/"^; 



M.3i[+ - -+] = ^'^^{FHml - ml)' + J^(F+ + 2mf)}V^ 



Mi3i[- - +-] = ^.^Smi[{ml - m\f + J^\F^ + 2mf (m? - m^)2}V4. 
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2mlmj + m^(-2mf - 3m^i^i + 6m^A)] + [3rr?4( + mf ) + m^(-3J^A + 3mf A - 

28mf J'^-SOml^mf J2-12mfmf)+32m|°i^2^+6mfmfA+mf(-28m2A+60m^)+2m^m^Ai^g- 
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Mi2i[ ] = -^--^{3m8mfF- + 2m\v4ml[{-2mf - A + m\)J'' - QF'] + 

6mfmfmj[-3F- + {-2Hf- \ + m\),P] - 2mjmfm'l[6m\mf + 32mf - ISmfm'^.Hf + 
Amf - 14m^A - Smfm'^A - 2m^m^ + m\{-2mfHl'^ - Sm^H^ - 6mfA)] + m^^[3m\{J^ + 
mf) + m^(3j8A - 3mf A - 28mp'^ - SOmfmf - 12m^m{) + 32mfH^ - Gmfrnf A + 
mf (28m2A + 60mf) + 2mfmf Aifg _ 2m^mfif|^]}^/^ 

M3i[- + +-] = - + J'V^'; 

M3i[+ + ] = ^^-^{mfmfF++2mfmfm2[j2(2m2-A-m2)-2F+]+2m|mfm^[3F+ + 

J^(3m^ + 2A — 2FPi^\ — 2rnimfm^[2mf{m\ + mg) +m^A — 2mf{A — 2m'l) +m? (-3m^A + 
2mf ) + m\{m^Hf - 2mfH^ + 2mfA)] + m^[(A2 - m^A) + mf (m^ + m^A - Bm^m^) + 
2m\mfmiHf - 2mfm2AiJ| + 4mfmfii"| - 2mfm^,AHf + 2m|mf]}^/^; 

M3i[ ] = ^^^{mfmfF-+2mfmfm2[J2(2m2+A-m^)-2F-]+2m|mfmf [3F- + 

J2(3m^ - 3A - 2i?2 J] - 2mjmfm^^[2mj{-m\ + ) +m^A - 2mf (A - 2m2) +mf(3mf A + 
2m^) + m\{-mlHl + 2m|i7| + 2m|A)] + [(A^ + m\A)J^ + mf (m^ - m^A - 9>m\m1) + 
2m\mfmjHf - 2mfmlAHf + Am'^^mfHl + 2mfm^^AHl + 2m^mf]}^/^; 

Mn[- + +-] = ^^(F+)V4^(^2_^.)2 + j4. 
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In the event that decays proceed through off-shell sleptons, the endpoint expres- 
sions simplify dramatically, e.g. for M^i, M 2121, M 121, M^i, and Mu, either [-|- -| ] 

or [ (alternatively [ ] or [-|- -|- respectively) are maximal. For 
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It should be easy to see the difference between this pattern of endpoints and those 
for the on-shell case. 



